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ABSTRACT 
The cross sections for the elastic scattering of electrons from 
atomic oxygen are calculated in the energy range between O and 13.6 
electron volts. Polarization of the atom is included by using a term 
in the potential function which has the correct asymptotic behavior and 
which vanishes at the origin. Previous theoretical and experimental 


results are discussed. 
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1. Siiteroduerron 

A knowledge of the elastic scattering cross sections for electron- 
oxygen collisions is of practical importance in the study of the pro- 
perties of air at high temperatures. The experimental determination of 
low-energy cross sections is especially difficult in the case of atomic 
oxygen because it recombines easily to form molecular oxygen. Only re- 
cently have experiments been done to measure this cross section in the 
energy range below 15 electron volts. [11,15] Several previous theo- 
retical treatments have been made. 11,3),7 OIDs | The results are 
divergent, and none agrees satisfactorily with the experiments. 

The elastic scattering of low-energy electrons by neutral oxygen 
is treated here using the partial-wave method of Faxen and Holtzmark.[ 13] 
This method, which is outlined in the next section, allows one to de- 
termine the cross section if the potential field in which the electron 
finds itself is known. The first approximation to the scattering treats 
the atom as a positive charge at r = O surrounded by a spherically sym- 
metric negative charge distribution. The charge density of this elec- 
tronic distribution is calculated from the Hartree wave functions for 
oxygen. 

This first approximation is inadequate because it ignores the dis- 
tortion of the atomic electron distribution by the Coulomb field of the 
incoming electron. This polarization effect is included in this calcu- 
lation by representing it with a polarization potential in addition to 
the undistorted potential. The polarization potential can be thought of 


as the term which represents the difference between the true interaction 


potential and the undistorted potential. The form it assumes for large 
Separations is known; for small separations the form is uncertain and 
will be discussed further in Section 3. 

When the incident electron is within the electron cloud of the 
atom there is an additional repulsive force due to the Pauli principle. 
Hammerling, Shine, and Kivel have treated this exchange effect by the 
inclusion of a third potential to approximate the change in the electron's 
energy due to the deficiency of electrons with like spin near it. [7] 
Such an approximation is similar to Slater's treatment of exchange in 
atomic wave function calculations. The result of their cross section 
calculation which includes this approximation shows that the exchange 
effect makes a negligible contribution to the cross section in the energy 
range considered here. Furthermore our polarization potential will have 
a parameter which is adjusted to fit experimental data. The value of 
this parameter is affected by the exchange contribution, so no specific 
account is taken of the exchange potential in this calculation. 

The actual evaluation of the cross sections using the phase shift 
method involves using a computer to integrate Schroedinger's equation 
numerically. The details of the procedure are given in Section 4. 

The following atomic units are uSed: 

unit of mass =m = rest mass of the electron 
unit of charge = e = charge of the proton in Gaussian units 


ao = ee ee of the first Bohr orbit 


unit of length 
, Z ; 
Consistent with these, the unit of energy is e [ass which is twice the 
ionization energy of hydrogen and the unit of angular momentum is h/2m. 
It should be noted that one atomic energy unit equals about 27.2 electron 
volts; this is twice the value of another commonly used unit, the rydberg, 


Wipen 1s 13.6 ev. fee 


2. Review of the partial wave method for computing the cross section. 
The scattering of two bodies when considered in the center of mass 
system is equivalent to the scattering of a single particle with the re- 
duced mass from a potential, V(r). For the electron-atom problem the 
reduced mass is very nearly equal to the mass of the electron. In this 
case the center of mass frame is identical to the lab frame. 
Since the potential function, V(r), is assumed to be spherically 


symmetric, the electron's wave function may be separated and written as 


¥ (r,9) Ry (r) Py (cos Q). 


In atomic units Ry (r), the radial wave function, satisfies the equation 


dR, (r) 
1 _d Mii) ; — £41) _ 
2 dr | wae + [2E - 2V(r) 2 ie) = 025) (1) 





If we assume that there is a distance, r = a, beyond which V(r) is zero, 


equation (1) may be solved exactly in the region r > a, giving 
neCL) = Aeeiees 66 jy, (kr) ™ simnee.n, (ke | (2) 
£ £ Loh Lh : 


where Ay and wr) are constants of integration, Jy and n, are the spherical 
Bessel and Neumann functions, respectively, and k = V2E. The asymptotic 
form of equation (2) is 


A 


> £L ; 
R(x) ~ ar [sin (kr = 4 in + 54) ip 
so that 
> A 
L 
¥ (r,0) ~ = P, (cos 0) — [sin (kr - 3} £0 +5,) ]. (3) 
2 kr L 


The physical meaning of v) may be seen by realizing that if there had 


been no potential anywhere, the asymptotic form of Ry (r) would be 


A 
ze Sin (kr - 4% £11). Hence ae the phase shift angle, is the quantity 


kr 
which describes the scattering. 

In a symmetrical scattering situation we expect an asymptotic wave 
function to be in the form of a plane wave superimposed on a spherical 
scattered wave. If the plane wave propagation is in the z direction, we 
express the asymptotic form of the wave function as 


ikz 


Yre +4 £(@) oS 


Since the asymptotic form of eee is 
kz « ,d (kr-} 2 
. ’ 
giz) i841) Ba(cos 0) SmeGereeetD) 
2=0 £ kr 


the asymptotic scattered wave function may be written as 


co 


e o~ & e 
LPO os Ge SS ee ae (4) 
2=0 L kr 


Equating (3) and (4), writing the sine function in exponential form, 


and solving for f£(@), we get 


eee OS cite 4-1) BGes 0) 
£=0 A 


The differential scattering cross section is 


aa) Z 
sem | £¢@)) | ‘ 


Integrating this over all solid angles, we get the total cross section 


oni Ee ty) sin’6, . (5) 


The value of 6, will be very small for £ > ka, so that for the range 


L 
of energies considered here values of £ greater than 2 or 3 need not be 
considered in calculating the cross section from equation 5. 

Hence in order to find 0, we need the values of ay . These-ane 

F < 

found by solving the wave equation (1) numerically to evaluate Ry (r) 

: : ; > , : ; 
and its derivative at r =a. Ry (r) and its derivative are evaluated 
at r =a from equation (2). These two wave functions are matched at 
the distance a by equating the logarithmic derivatives. This gives 
the values of the phase shift angles, 5 When all significant phase 
shift angles have been calculated, equation (5) is used to evaluate the 
cross section at the given value of energy. This process is repeated 


over the energy range of interest to find the spectrum of cross sections 


versus energy. 


3 Potential function for electron-oxygen scattering. 

A calculation of the electron-oxygen elastic scattering cross 
section is routine, although tedious, if V(r) is known and a computer 
is available. In this calculation we assume V(r) is the sum of the 
undistorted atomic potential, VX), and the "polarization" potential, 

V fe)» which attempts to represent the difference between the true po- 
tential and the undistorted potential. 

Ve) may be found in a direct electrostatic manner if the electron 
charge distribution is known. [14] The charge distribution is calculated 
from the Hartree wave functions of the oxygen electrons. Since the 
Hartree approximation treats each atomic electron as moving in a spheri- 
cally symmetric potential, the atomic electron wave functions can be 
written as 

Unt = Rng (t') ¥g,(0' 8"), 
where Yance 29") are the normalized spherical harmonics. The values of 
Rig it? have been tabulated by Hartree. [8] In computing the electron 
volume charge density one sums uu for each electron. For filled n,& 
shells this sum will be spherically symmetric; the only asymmetry will 
be contributed by the outer shell which is not filled. It is within the 
Hartree approximation to neglect this asymmetry. Hence the charge density 


of electrons is given by 


Zs 
(\ Se ' 
O(r‘) =2 Nu R g(t ) 
n,& 


where Ni represents the number of electrons in the n,£ shell. Only the 
3p ground state of oxygen is considered here, for which Nio = 2, Noo = 2, 


and Noy = 4, 


=10- 


V(X) consists of a negative term due to the nucleus and a posi- 


tive term due to the electron cloud. It may be written as 


Di DT 


2 2 2 
— P2°Re (ry + ZR, Cr ae Ar +) 
Vr) ee + , ee me Q@' dO@'d@'dr'. 
Oe Ue | rc | 


Letting riR (x ) = Pg tt Dis 


roe) 2 
[P (c')+ PS Fait ee dg) 
v(t) = Beonf aa FaalOe Me I oe 
0 r- r' | 


The integration can be divided into the contribution of the electrons 


between r' = 0 and r' = r, and the contribution of those from r' =r 
to r' =, Then Via may be written as 
Lente’ ) RIES) Re Cel 
_ 8 10 20 sig ; 
Vr) =- = + 8t SS ee (6) 
H r r 
0 


+m f Ma [Pi (r')+ "06" yee te) Tl) 
dr 


2 


P 1 P , : 
Ns ng St ) is called the radial charge density. It 


The function &% 
n,& 


represents the probability per unit length of finding an electron in the 
spherical shell between r' and r' + dr'. Figure 1 is a graph of this 
function for oxygen. 

The treatment of the polarization potential is less exact. The pro- 
cess of treating the atomic potential as the superposition of an undis- 


torted potential and a polarized dipole potential is less than rigorous. 


- we 


However, even if we accept this approach, the form of the polarization 
potential, V(r), is not Certain. 

lf the incoming electron is at large r, its electrostatic field at 
the atom is nearly homogeneous and has a value of “ie, This field in- 


duces a dipole moment which is directly proportional to the field, 


Lb = ye where @ is the polarizability. This dipole moment then exerts 


an attractive force back on the electron with a value of £ = Bye = 2etae 


So, for large r, 


lo ©) 
5 a oy 
V(r) ~-f 2¢/r" dr = - ae 
r 


This form of Vp (x) for large r must be modified so that the function 


does not blow up at the origin. At the origin there should be no polari- 


zation force since, in this limit, spherical symmetry prevails and no 
multipole moments can be induced. For this reason the proper polari- 


zation potential should have the property that 


Limit d _ 
ae ae 


When the scattering electron comes "inside" the atom, the atomic elec- 
trons still tend to hide behind the nucleus. However in this case the 
asymmetry gives rise to a repulsive force. In this classical picture, 
the polarization potential should also have zero slope at a value of r 
which is characteristic of the size of the atom. 

A function which satisfies these conditions and which will be used 


in this calculation i8 


e4 


Vp (r) = ee 1 - exp [- im, ¥ ey (7) 


m= 


This function, with q equal either five or eight, was first used by 
Biermann to account for polarization in the calculation of oscillator 
line strengths. [2] Garrett and Mann have used this polarization potential, 
nr q = 8 in calculations of electron scattering from alkali atoms. [5,6] 
They obtained good agreement with experimental scattering data by treating 
@ and to as parameters. Graphs of equation 7 are given in Figures 3 and 4. 

The exponential in equation (7) may be expanded so that for small 
r we have 

7 0 ald 

From this we see that the requirement that the polarization force vanish 
at the origin restricts q to be greater than 5. In this calculation we 
use values of 5, 6, and 8 for q. The result with q = 5 represents a 
limit since any value of q larger than five satisfies the condition 
mentioned above. 

The value of ty is taken to be the expectation value of the radius 
of the valance electron shell. Its value is calculated from the Hartree 


wave function for the outer shell as follows 
co 


ry |. rag Pee dnis 
This yields ae 1.2 atomic units. 

No direct experimental values of the polarizability are available. 
Several calculations based on indirect experiments have been made. Klein 
and Brueckner found @ to be 5.589 by adjusting it to give the correct 
binding energy for 0. [9] Using a similar procedure, Cooper and Martin 
calculate a value of 5.499. [3] Both these methods do not attempt to 


provide the true polarizability of the atom, but yield a number which 


aie A 


includes exchange and other effects. Dalgarno and Parkinson found a 
value of 6.005 based on experimental determinations of the dielectric 
constant. [4] In our treatment to is fixed and @ is taken to be either 
5.5 or 6.0, and for both of these values the cross section is computed 


for q = 5, 6, and 8. 
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4. Numerical Methods. 

Before the wave equation can be integrated it is necessary to know 
V4). This function was calculated from equation (6) using a trapizodial 
summation technique. It was found that the resulting function could be 


fitted closely by 








ia “Aor -\ x - Yr 
V(r) =- 8 Ly = + ¥y = +v., e . aC > ] 
H Loar 2 r 3 4 : 
where Y) = 1.2806 Your -0.2806 
dy = 2.2376 ro = 18.263 
V3 = 3.0715 We = -1.9710 
de = 6.803 di, = 11.548. 


This functional form was much easier to enter into the computer than 
tabulated data points, and the fitted function matches the computed ~™. 
function the fifth significant figure. 

To solve equation (1) the substitution Y, (x) =r Ry () was made, 


yielding the equation 


{ 


2 

d Y,(r) 

- + = (pene yy 2V(r) + 2, (2) - 2B } ¥, (x). (8) 
dr iG 


Equation (8) was integrated using a Runge-Kutta approximation. In order 
to use such a method it is necessary to know an initial value of both Y 
and dY/dr at some starting point. These values were found by assuming 


asolution 


= 152 


Since this solution will be used at a small value of r, a series approxi- 
mation for Vo) and Vp (x) may be used in equation (8). Neglecting terms 


: 2 
in r or smaller, we have 


pie pew bo/r 4 b, + bor, 


where 


om 
' 


VM Yoh5 
by = S| Suet n>) T Y3 Aamuaye) 








Altsomforesiall reev as —~ Only for q = 5 will this make a con- 


po 


tribution to the total potential in the above approximation, changing the 





value of b, given above to b, + ; » Equation (8) for small r, becomes 


c 27 
fe) 





2 a°y Z 3 Z 
r +1 - £(4 +1) + 2(b) 2 ae bor + bot + 2E r) eo = 0. (8') 


Substituting the series solution for Y. into equation (8), we get the 
indical equation a@(@ - 1) = £(4£ +1). The only acceptable root is 


w= +1 since Y(0) = 0. The recursion relation is 


of 
. aa Ly ey arial? OE aes Soc pag 
= mar ee + 1) 
dY (r) 
£+1 8 Late 2 s L, Sita) of + 1 
So Yj (c) wr > ee, and ~w (£+1)r° + tet ; 
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One could, of course, evaluate these equations at r = O and start the 
integration from there, however it was found that the accuracy of the 
Runge-Kutta integration approximation is very sensitive to the step 
size chosen near the origin. For this reason we evaluate Y. and d¥ /dr 
at r = .Ql and use these values to start the iteration procedure. The 
step size was determined by choosing an arbitrary value to do a trial 
calculation, reducing the step size and repeating the calculation until 
there was no Significant difference between two successive results. A 
step size of .002 was used. Although this value is smaller than needed 
for large r, no attempt was fade to change the step size as the iteration 
progressed. 

Equation (8) was integrated to the final point r = a which was 


determined by the condition 


Vp (a) + Vi (a) yi 
cy ore 10 


For r > a we have a region of almost zero potential, so from equation 


(2) the solution of the wave equation is 
> — e s 
Y,(r) = Ayr [ .cos 8 jy (kr) - sin 6 yn, (kr) ). 


The derivative is 


> 

dY) (x) k 

— Ayr aT f Cig. Cee) - (4 + 1) Jg4, (kr) | cos 6, : 
Ye) 





-[én,_, (kr) - @ + 1)n,, (ker) sin 5 + 


‘ < <a. 
Letting Y, (a) and Y) (a) denote the computed values of me (r) and 


ts 


== 


av) (x) 


ce at r =a, we equate logarithmetic derivatives at r = ato get 


<! 

Yo (a) k 
a eS a a SS [4 j,_, (ka)- G41) j, | (ka) Jeos 7) 
Y, (a) (2£+1)[ cos 6 pj, (ka)-sin 6 yn, (ka) | 


Ling (ba) Wag (a) to 





Solving for o gives 
k Y; (a) 
a | Ziel Hy. 1(ta) ~ GHD yyy ha) ] “2 ig (ead} + 3, (4a) 
bY = arctan Ce (a 
k ¥, (a) 
«(opal An, _, (ka) . (241), (ka) | - a ' ny (ra) + n, (ka) 
Y, (a 
L 


From equation (9), for a given £ and E, a principle value of 6, is found. 
For any E, this process is repeated for £ = 0, 1, 2, and 3. Then the 
value of O© is found from equation (5). As expected the contribution to 
the cross section from the £ = 3 term is negligible in the energy range 
considered. 

The entire calculation was done for energies between 0 and .50 in 
steps of .05 for q = 5, 6, and 8; and for @ = 0, 5.5 and 6.0. All con- 
putations were done on a C.D.C. 1604 computer using Fortran-60. The com- 
putation time for each run (a given q and @) was about four hours. This 


large amount of computer time could be reduced by increasing the step 


size at large values of r. 


See 


5. Previous calculations and experimental data. 

The theoretical treatments of the low-energy scattering problem 
may be divided into two groups. One approach is similar to the one 
used here, but different forms of the polarization potential are used. 
The other treatment is much more elaborate. It attempts to solve the 
Hartree-Fock equations of the system of scattered electron plus atom. 
Figure 7 shows some results of previous calculations. 

Robinson uses no polarization potential. [12] His potential function 
is calculated from Hemeiie anatase functions, and the cross sections 
calculated by the partial-wave method for 2 = 0and1. He also evaluated 
the cross sections using the W.K.B. approximation and showed that it leads 
to very large cross sections at low energies. 

Klein and Brueckner also treat the problem using the partial-wave 
method. [9] They account for polarization with the potential 


my 


2h + oe 


V(x) = ; 
This potential function has been used often in the treatment of polari- 
zation in atomic wave function calculations. It gives a zero polarization 
force at the origin, but differs from the function used here (equation 7) 
in that it produces an attractive force for all r. As mentioned before, 
they determined the polarizability from experimental data on the binding 
energy of O . Their curve in Figure 7 includes’ s and p waves (4 = O and 
1). Cooper and Martin [3] found a mistake in the paper by Klein and 
Brueckner. They use the same procedure and arrive at cross sections which 


show a decrease at low energies. 
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Hammerling, Shine, and Kivel treat polarization in the same way as 
Klein and Brueckner. They take @ = 5.4 and r, as an adjustable para- 
meter. [7] They also include a term in the total potential to represent 


3/2 


the exchange energy. This term varies as l/r and is derived from 
the Slater approximation. They show results obtained by neglecting both 
polarization and exchange, and neglecting polarization only. The cross 
sections obtained by using just the Hartree potential are different from 
any obtained previously. The cross sections with exchange and polari- 
zation are plotted in Figure 7. It is not clear why these should be so 
different from the other calculations since they reported using the same 
method. 

Temkin includes polarization effects within the Hartree-Fock for- 
malism. [16] He associates with each atomic electron's wave function 
a polarized part which depends on the coordinate of the scattered elec- 
tron. The wave function of the scattered electron is made antisymmetric 
with respect to the total atomic wave function. This gives rise to the 
kind of polarization potential considered here, called direct polari- 
zation, and another term, called exchange polarization, which contributes 
about 15% to the cross section. His results for s-wave scattering are 
shown in Figure 7. Bauer and Browne have modified Temkin's method. 1) 
They made more drastic simplifications in the Hartree-Fock equations and 
are able to calculate s, p, and d-wave scattering. Their cross sections 
are small at low energies and fairly flat at higher energies. 

The most accurate experimental data is that of Neynaber, et al. [11] 
They have measured the cross sections in the energy range from about .1 


to .45 (2.3 to 11.6 eV). They find the cross section to be almost constant 
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at 19.5 + 1.3 throughout the energy range. This represents a band in which 
half of their 64 experimental points fall. Sunshine, Aubrey, and Bederson 
measured the cross sections between about .03 and .45. [15] There is 

more scatter in their data, but there is general agreement between the 

two experiments. Lin and Kivel obtained a low energy point at about .02 
energy units using a shock tube. [10] If this point is accurate it in- 
dicates that the cross section decreases sharply at low energies. Al- 
though these experiments are extremely difficult and have low precision, 
they do indicate the magnitude of the cross section and show its general 


dependence on energy. 
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Gs Results. 

The numerical results of our calculation are given in Tables ous 
and 3 and summarized in Figures 5 and 6. The curve labeled a = O repre- 
sents the cross sections obtained by ignoring the polarization potential. 
This curve corresponds almost exactly to that of Robinson. This indicates 
that the differences between Hartree - Fock wave functions used by him 
and Hartree wave functions used here have an insignificant effect on the 
cross section. 

The introduction of the polarization potential increases the phase 
shift angles but decreases the cross section. The smaller values of q 
produce potentials which give rise to smaller cross sections, and also 
straighten the curve out. The remarkable feature indicated by Neynaber's 
experimental results is the almost constant value of the cross section 
between .1 and .45 energy units. [11] Our results give a similar constant 
value in this energy range. Further, the measurement of Lin and Kivel 
gives evidence that the cross section has a very small value at the 
lowest energies. [10] This is also a property of our results. 

Of all the theoretical calculations that have been discussed, that 
of Bauer and Brown is the most pure (not dependent of experimental para- 
meters) which extends of the entire energy range. [1] Our calculation 
with @ = 6 and q = 6 almost duplicates their results. Our results 


are also consistent with those of Temkin at the lowest energies. [16] 
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7. Conclusions and acknowledgements. 

No attempt has been made to justify the fundamental validity of 
the form of the polarization potential except in a general way. Experi- 
ments indicate that consideration of polarization in any scattering cal- 
culation is important. Although the "correct" way to treat the effect 
has been outlined by Temkin, it is desirable to have a more tractable 
semi-phenomenological method. Our results indicate that the potential 
function used here (equations 6 and 7) with @ = 6 and q = 5 gives cross 
sections which have the energy dependence indicated by experiment and 
by the more exact calculations. Our results seem much more realistic 
than previous treatments using other forms of the polarization potential. 

The author would like to thank Professor William B. Zeleny for 
suggesting this subject as a thesis topic. The help provided by the 
staff of the Computer Facility of the Naval Postgraduate School is 


greatly appreciated. 
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TABLE 1 ALPHA = 5.5 


Q = 5 
5 DEL TAQ DELTA1 DELTA2 DELTAS SIGMA 
~0050 1.9770 .5018 ~0000 -0000 6.6678 
~0500 1.8872 29104 ~0055 -0017 15.7451 
»1000 1.8269 2593558 .0117 »0035 19.4121 
.1500 1.7813. 2561 .0188 -0053 21.5178 
.2000 1.7440 7) Sey .0267 -0072 22.7739 
.2500 Lee 059955 ~0354 -0091 23.5334 
~ 5000 1.6843) -6091 -0448 ~0112 24.0109 
»3500 1.6594 -6230 ~0550 »01352 24.3357 
~ 4000 1.6371 »-6356 -0658 -0154 24.5858 
»4500 1.6167 26471 ©0771 -0177 § 24.8064 
»5000 1.5979 clei! , .0888 -0200 25.0224 

Q = 
E DEL TAO DELTA1I DELTA2 DELTA3 SIGMA 
»0050 1.9709 .5006 ~0000 - 0000 10.4912 
.0500 1.8719 9274 ~0055 -0017 22.2738 
~1000 1.8078 -9611 »0117 »0035 27.5908 
©1500 1.7598 ~5889 .0188 -0053 30.0143 
»2000 1.7209 6121 .0266 -0072 30.9885 
.2500 1.6878 .6318 00353 ~0091 32.2518 
~ 3000 1.6591 -6488 00447 -0112 31.1230 
~ 3500 1.63355 ~6638 »0547 -01338 30.8598 
~ 4000 1.6105 ©6771 .0653 »0154 30.5507 
~ 4300 1.5897 -6890 ~0765 -0177 30.2516 
~5000 1.5706 -6998 .0880 20201 29.9895 

Q = 
= DEL TAO DELTA DEL TA2 DELTA3 SIGMA 
~0050 1.9644 .59028 .0000 -0000 15.9588 
~0500 1.8554 »59582 ~0055 -0017 36.8699 
~1000 1.7871 »6009 ~0117 -0035 42.9553 
»1500 1.73566 -6318 ~0187 -0053 43.8549 
~2000 1.6960 »-6562 ~0265 -0072 43.0616 
»2900 1.6617 ~6762 »0351 -0091 41.7041 
» 3000 1.6320 -6934 0444 ~0112 40.2218 
~ 5300 1.6057 ~7083 .0543 -0133 38.7979 
~ 4000 1.5822 ~7214 -0648 -0154 37.5062 
~ 4500 1.5610 -7332 .0757 »0177 36.3708 
»5000 1.5416 .7438 ~0869 -0201 $5.3891 


aoe 
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.0050 
»0500 
»1000 
~1500 
~2000 
»2300 
» 3000 
» 3500 
»4000 
©4500 
»5000 


TABLE 2 

Q = 
DELIAO DELTA1 
1.9825 22023 
1.8977 .2040 
1.8388 2234 
1.7939 ~5434 
1.7569 .5619 
he 7295 72786 
1.6976 9937 
1.6729 ~-6074 
1.65905 -6198 
1.6301 -6312 
1.6114 .6416 

Q = 
DELTAO DELTA1I 
1.9756 .5014 
1.8810 ~2186 
1.8161 22483 
1.7707 5745 
1.7320 -9970 
1.6992 ~6165 
1.9670 6336 
1.6450 6487 
1.6220 -6622 
1.6012 6744 
1.5820 6854 

Q = 
DELTAO DELTA1 
1.9665 -5006 
1.86382 -5426 
1.7959 -9830 
1.7458 -6140 
1.7054 -6390 
1.6712 ~6599 
1.6415 ~6778 
1.6152 -6934 
1.5917 ~7072 
1.5703 » 7196 
1.5509 ~7307 
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ALPHA 


= 6.0 


DELTA2 


~0000 
-0060 
-0129 
.0207 
.0294 
»0390 
-0496 
-0609 
.0729 
-0856 
.0967 


DELTA2 


~0000 
-0060 
»0128 
» 0206 
-0293 
.0389 
.0493 
.0605 
.0724 
» 0848 
.0976 


DELTA2 


-0000 
.0060 
-0128 
»0206 
»0292 
-03587 
-0490 
- 0600 
~0717 
»08358 
-0963 


DELTAS 


-0000 
~0019 
-0038 
-0058 
»0079 
-0100 
»0122 
-0145 
-0184 
-0193 
-0218 


DELTAS 


-0000 
-0019 
-00358 
-0058 
-0079 
-0100 
.0122 
©0145 
-0168 
©0193 
-0219 


DELTAS 


-0000 
-0019 
-00358 
-0058 
-0079 
-0100 
-0122 
»0145 
-0168 
0193 
0219 


SIGMA 


4.0100 
12.8560 
16.3795 
18.5432 
20.0096 
21.0657 
21.8826 
22.5658 
23.1803 
23.7635 
24.3541 


SIGMA 

3 

7.4262 
18.2930 
23.2639 
25.9655 
27.4278 
28.2089 
28.6278 
28.8643 
29.0231 
29.1604 
29.3036 


SIGMA 


12.2868 
28.7581 
35.5747 
37.9080 
38.3704 
38.0560 
37.3904 
36.6634 
35.9643 
35.3414 
34.8102 
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E 


.0050 
»0500 
»1000 
»1500 
-2000 
-2900 
$000 
» $500 
» 4000 
» 4500 
»5000 


TABLE 3 


DELTAO 


1.9154 
1.7502 
O67 2 
1.6121 
eo 205 
1.5370 
1.5091 
1.4850 
1.4638 
1.4450 
1.4280 


DELTA 


.2004 
me 157 
29498 
-6000 
-6664 
~7415 
»8139 
~8750 
92238 
~9575 
»9836 


ALPHA 


= 0.0 


DELTA2 


.0000 
-0000 
-0004 
-0011 
»0020 
»0032 
»0047 
-0064 
-0084 
»0105 
.0128 


DELTAS 


-0000 
~0000 
-0000 
-0000 
-0000 
»0002 
-0004 
~0005 
»0008 
-0011 
-0014 


SIGMA 


86.6421 
63.6804 
51.5073 
48.9037 
53.3931 
60.4892 
64.6529 
635.9188 
59.9304 
54.7727 
49.6356 
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